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ABSTRACT 


The differential cross sections for elastic scattering of 180 MeV protons from lithium, beryllium, 
carbon, aluminium, calcium, iron, indium and gold have been measured up to large scattering 
angles. Particular care was taken to obtain good energy resolution. The results have been com- 
pared with the calculated scattering from an optical model potential with different radial exten- 
sions of the real and imaginary terms in both the central and the spin-orbit parts of the potential. 
Qualitative agreement was found. 
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1. Introduction 


The differential cross sections for elastic scattering of high energy protons from 
nuclei have been measured in many experiments ever since the first beams of such 
protons were available from accelerators. On most occasions, however, the measure- 
ments have been confined to small angles and in some cases the energy resolution has 
been insufficient to make sure that the cross sections determined were unaffected by 
inelastic scattering. 

Many analyses of high energy proton scattering data have been made using optical 
model potentials. Several reviews have been written [1, 2, 3], where references are 
given both to theoretical and experimental work. Usually the real and the imaginary 
central terms of the potential have been assumed to have the same radial dependence 
and the spin-orbit terms to be proportional to the gradient of the central parts. 
However, such a potential is insufficient to explain the scattering from com- 
paratively light nuclei [4, 5, 6]. A more general potential with different form factors 
greatly increases the number of adjustable parameters in a phenomenological fitting 
procedure and consequently requires more experimental results for a successful 
analysis; large angle data should be particularly useful. 

Under certain assumptions it has been possible to reduce the actual many-body 
problem of the scattering of a high energy proton from a complex nucleus to the 
problem of the scattering of a proton from a generalized two-body potential, which 
however is non-local for finite nuclei (see for instance [7]). In a local approximation 
the optical model potential is obtained from the nucleon-nucleon scattering matrix 
and the nuclear density. If the angular dependence of this matrix 1s taken into 
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arrangement, 01 only very few smali-angle cross sections were d 
The second part fad the paper consists of an optical model analys 


cross sections together with reaction cross sections reported in ‘the ace 
paper [8] and previously published polarizations [6, 9-13]. Two types of fits | 
been tried, one purely phenomenological and the other using the Gammel-’ 
nucleon-nucleon potential and nuclear densities derived from electron HORUS 


[14]. . 
2. Experimental 


2.1. General considerations 


The main experimental problem in determining the differential cross sections for 
elastic scattering from nuclei at large angles is to obtain sufficiently good energy 


resolution. In the present experiment essentially the same method was used as tha 
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Fig. 1. Plan view of the experimental arrangements. 
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by Lyx n and Maris [15] to measure 
e. the 


differential cross sections for inelastic 
the separation of elastically and inelastically scattered protons was 
red by using a magnetic spectrometer which could be rotated about the scat- 
er he experimental lay-out is shown in Fig. 1. see teh 
The incoming beam was the external, unpolarized beam of 187 MeV protons from 
Uppsala synchrocyclotron. The energy spread was not known but was definitely 
s than 1 MeV full width at half height. 

The protons that emerged from the spectrometer were detected by a counter tele- 
ype of three scintillation counters in coincidence, the width of the defining crystal 
being so small that it did not appreciably reduce the resolution. The incident beam 
was monitored by another telescope looking at the end window of the proton tube. 
; The area under a spectrum peak corresponding to elastically scattered protons is 
proportional to the elastic differential cross section. The proportionality factor is very 
geently independent of the peak width if the momentum or the energy of the protons 
is taken as abscissa, but it is still a function of the magnetic field and the proton 
energy. The determination of this factor is considered in section 2.5. and was based 
on values from other experiments for the free p—p scattering cross section. 


it 
2.2. The magnetic spectrometer 


The main properties of the magnetic spectrometer have been described by Tyrén 
_and Maris [15]. Some changes of its field shape were made, however, to improve the 
accuracy of the instrument. 
_ The magnetic field was not quite homogeneous in the whole region where protons 
- passed due to the small width of the pole tips. (The gap was 75 mm wide and 16 mm 
high.) This would probably reduce the energy resolution and, what would be more 
‘ serious, make the transmission of the protons dependent on the proton path inside 
the spectrometer. It was not practical to use such a small height of the entrance slit 
that all protons were prevented from hitting the pole faces and defocusing (and fo- 
cusing) forces due to magnetic field gradients would thus strongly affect the trans- 
‘mission. Therefore, the field shape was altered at two places along the proton path by 
shimming in such a way that the direction of the gradients were reversed. The effect 
of these shims was then to make every proton pass through successive regions of 
vertically focusing and vertically defocusing forces, the net effect of which was 
approximately zero. The energy resolution should also be improved by this shimming 
because the average value of the magnetic field experienced by a proton was made 
nearly independent of the path. 

The second change consisted of a careful shaping of the fringing field both at the 
entrance and at the exit to make it satisfy the theoretical requirements more ri- 
gorously. 

As a result of these changes the following improvements were found experimentally: 

(1) The transmission of protons entering at different points of the opening slit, 
which was 8 mm high and 12 mm wide, was constant within a statistical uncertainty 
of a few per cent. 

(2) Within 5 % the transmission of protons was found to be independent of changes 
in the beam position at the scatterer of + 1.5 mm vertically and + 10 mm. horizon- 

tally. This was quite sufficient to make sure that the errors that could arise from small 
beam displacements or from the use of thick targets at large scattering angles were 
negligible. 
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Tt was found that the extensive fringing field gave rise to the beam | 
sition at the scatterer, a shift that varied not only with the magnetic field but 
with the-angular position of the spectrometer. By small adjustment magnets 
Fig. 1) about two meters in front of the scatterer the beam was always brought 
the same position, and corrections were applied to the zero degree correction an 
Immediately before and after the measurement of a peak, the beam position at the 
scatterer was determined by exposing a photographic paper. It was thus ensured 
that no changes had occurred either in the vertical direction, which might have al- 
tered the transmission, or in the horizontal direction which would have deformed the 
peak shape. ¢ 
--~ Jn a scattering experiment the scatterer is usually set at such an angle that the 
normal to the target surface bisects the scattering angle. All protons then have equal — 
path lengths in the target and they emerge with approximately the same energy. — 
In this experiment the target angle was chosen more elaborately to obtain as small 
peak widths as possible. The following sources of energy spread were considered and 
eliminated to first order by a proper choice of target angle: (a) the differences be- 
tween the recoil energy loss in the beginning and the end of the target, (b) the non- | 
linear variation of ionization energy loss with proton energy, and (c) the width of the © 
beam as seen from the spectrometer. , 


2.3. Electronics 


i aT 


Due to the fact that the duty cycle of the external beam from the cyclotron was 
very small (about 1/400) there was a risk of counting losses even at moderate count- 
ing rates. Moreover, the duty cycle was not very constant in time which would have 
made corrections for the counting losses due to dead time very uncertain, so that fast — 
counting systems were needed. 

Two scale-of-four units were built with fast transistors (2N 501) and used in the 
monitor and the spectrometer detecting systems. The resolving time for two conse- 
cutive pulses was 25 nsec and consequently two protons appearing in two consecutive 
rf. bursts of the proton beam could be resolved. The pulses from the fast scalers 
were fed into conventional 1 Mc sealers. 

The systems of photomultipliers, coincidence units and scalers were tested with 
varying counting rates, and counting losses were found to be less than 1 % for count- 
ing rates up to 10% protons/sec average rate. The number of protons detected per 
second was always in the experiment kept less than 10%, so that errors due to dead 
time losses should be negligible. 


2.4. Energy calibration 


A calibration curve for transformation from spectrometer current to proton 
energy was obtained by measuring the current at the maxima of a series of peaks 
from p—p scattering at several angles and from elastic and inelastic scattering from 
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anticoupled coils in series with each other and with a sensitive fluxmeter. The 
ere placed in the fringing field of the spectrometer and one of them between 
o mumetal bars. Each coil was shunted by an external resistance and the actual 
g consisted of an adjustment of one of the shunt resistances to give zero deflec- 
n on the galvanometer when the two coils were flipped simultaneously. By a 
table choice of coil sizes the device was made extremely sensitive to changes in the 
agnetic field in the range necessary for the experiment. The resistance readings were 
ut on an absolute scale each time a peak was measured by comparing with current 
eadings and using the calibration curve mentioned above. Thus the errors from the 
uncertainty in the energy setting of the spectrometer is small in most measurements. 


2.5. Transmission and efficiency calibrations 


im To be able to convert the area under a spectrum peak into an absolute value for the 
_ differential cross section one should know the following parameters: 

_ (a) The transmission of protons by the spectrometer magnet. 

_ (6) The effective solid angle for the scattering. The fringing field of the magnet 
extended outside the entrance slit and consequently the effective solid angle was 
- dependent on the focusing properties of the field. 

_ (c) The detection efficiency of the spectrometer counter telescope. Between the 
‘ second and the third counter of this telescope was inserted a copper absorber to pre- 
vent particles of low energy from being detected. Such particles could be protons 
having passed through the spectrometer by scattering on collimator edges or pick-up 
deuterons of the correct momentum to pass freely through the magnet. The efficiency 
was thus less than one mostly due to losses of elastic protons in this copper absorber. 

(d) The equivalent width in proton energy (=the abscissa of the spectra) of the 
defining crystal of the telescope. 

_ (e) The conversion factor between monitor counts and number of protons hitting 
the target. 

The first four of these parameters are dependent on the magnetic field, i.e. on the 
energy of the elastically scattered protons. 

Each of these five parameters was not measured separately. By determining the 
area under the peak for free proton-proton scattering from polythene minus carbon 
‘at some angles and using the relatively well known results for the p—p differential 
cross section, a calibration curve could be constructed for the product of all the five 
parameters versus proton energy. One additional advantage of this method was that 
any error in the energy calibration curve was eliminated from the end results. The 
reduction in beam intensity due to inelastic collisions in the target were corrected 
for. The values for the p—p cross section, which were used, are given in Table 1 and 
were obtained by an interpolation between published data [16-19]. 

Since the main absolute error of the cross sections reported in this paper originates 
in the uncertainty of the p—p cross sections, a future improvement of the latter would 
also improve on the results given in Tables 2-9. 
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about 0.01 MeV. Although the operation of the spectromet 
egion of strongly saturated iron, the ordinary t instrument 


For this purpose a device was used mich do 
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A. JOHANSSON et al., Elastic Scattering of Protons from Nuclei 


Table 1. The values used for the p—p differential cross section at 180 MeV laboratory 


energy. 

Lab. angle CM angle CM cross section | Lab. cross section 
(degrees) (degrees) (mb/sterad) (mb/sterad) 
8.5 17.8 3.74 16.16 
10.0 20.9 3.73 16.02 
11.8 24.7 3.75 15.98 
13.3 27.8 3.76 15.89 
14.3 29.9 3.76 15.80 
17.0 35.5 My SasrTe 15.48 
19.4 40.5 3.72 15.07 
21.7 45.3 3.68 14.61 
23.7 49.4 3.64__ 14.18 
25.5 53.1 3.61 13.80 
27.1 56.4 3.57 13.40 


The long time constancy of the calibration curve was checked by frequent remea- 
surements of a few carbon cross sections. Except on one occasion when a short circuit 
had occurred in the spectrometer no significant changes were found. 


2.6. The angular resolution 


As the elastic differential cross section varies extremely rapidly with angle nearly 
everywhere, it is of great interest to know the angular resolution function of the 
scattering, i.e. the resolution including the angular spread of the incoming beam, 
the beam width at the scatterer, the width of the entrance opening, and the multiple 
scattering in the scatterer. Therefore, the spectrometer telescope was replaced by an 
ionization chamber sufficiently large to cover the whole emerging beam. The spectro- 
meter was swept through the zero degree position with a few different targets at 
the usual scattering centre. In this manner were obtained not only the shape of the 
resolution function and its dependence on the multiple Coulomb scattering in the 
targets, but also the true spectrometer zero degree position which differed about 
0.7° from the nominal setting due to the extensive fringing field of the magnet. 

The experimental points for the resolution function were found to lie on a Gaussian 
with a total width at half height of 0.8° when the effects of multiple scattering were 
subtracted. The angular resolutions (with multiple scattering effects included) are 
given in Tables 2-9. 

There was a finite possibility for protons to be scattered twice in the target through 
wide scattering angles. The angular distribution of this double scattering would be 
expected to follow that for single scattering except for the absence of sharp diffrac- 
tion patterns. The effect would thus be most serious at the diffraction minima. An 
estimate was made which showed that double scattering could have contributed at 
most a few per cent to the elastic cross section obtained with the thickest targets used. 


2.7. Determination of peak areas 


A method was developed which permitted a convenient evaluation of the area 
under a peak known at some discrete points from experiment and also made it pos- 
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to estimate the effect of inelastic scattering to low-lying levels. The forms of 
spectrum peaks from “monoenergetic” groups of protons, i.e. protons without inelastic 
mixtures, were deduced as a function of half width by studying the peak shapes in 
he scattering at small angles from a few different targets, in the scattering from car- 
mn. at several angles and in the free p—p scattering. The peak shapes were found to be 
ussian except for a very slight skewness. Only curves with such shapes were used 
trying to fit graphically the experimental points of a peak plotted on a logaritmic 
cale. The area was obtained by multiplying the product of peak height and half 
width of the best-fitting curve by a known constant. The uncertainty in the method 
was estimated to be less than one per cent. 

_ For most scatterers used the shape of the elastic peak may be deformed on the 
low-energy side by inelastic scattering to low-lying levels. The excitation energies 
of these levels were assumed to be known. In the fitting procedure peaks of the same 
‘shape and width were used at the proper energy positions and with such magnitudes 
that the sum of all peaks gave a good fit to the experimental points. Although the 
inelastic peaks found in this manner were often very uncertain and the resolving into 
different peaks ambiguous in many cases, the uncertainty introduced in the deter- 
mination of the elastic peak area was comparatively small. A few examples of the 

spectra and the curve fitting are shown in Figs. 2-5. 


2.8. Inelastic contributions 


In Li’ the subtraction of the contribution from the first level at 0.478 MeV was 
quite uncertain due to the proximity of the ground state, but it was relatively small 
-at all angles and consequently the elastic cross sections were not affected very much. 
As all targets were of natural isotopic composition the recoil energy loss provided a 

“means of separating the scattering from the different isotopes if the scattering angle 
was large enough. Thus the Li-results refer to scattering from the natural isotopic 
mixture of Li (92.5% Li’ and 7.5% Li®) up to angles less than 25° and to Li’ alone 
for larger angles. Fig. 2 shows the Li spectrum at 27.8°. The elastic peak from Li® is 
included assuming the same cross section as for Li’. From the sum curve it is seen 
that the Li® cross section probably should be reduced somewhat and that there is 
some indication of inelastic scattering to the 0.478 MeV level. 

For Be the peaks from elastically and inelastically scattered protons were clearly 
separated by the spectrometer. Fig. 2 shows the Be spectra at 12.8° and 40.3°. 
Besides the strongly excited 2.43 MeV level there is also a contribution from a level 
or from the three-body break-up of Be at 1.75 MeV. 

The first excited state in C¥? at 4.43 MeV is too far from the ground state to contri- 
bute to the elastic peaks. Natural carbon has a 1.1 % admixture of C8, The elastic 
peaks from C!2-and C18 were sufficiently separated from each other at angles larger 
than 55° to permit the cross section to be evaluated for pure Cl’. Inelastic scattering 
in C13 was not considered. 

Al?7 has several low-lying excited states, viz. at 0.842 MeV, 1.013 MeV, 2.213 MeV, 
2.732 MeV etc. The fit to the inelastic scattering was therefore rather ambiguous as 
ean be seen from Fig. 3, but the elastic scattering was nevertheless quite well deter- 
mined. 

In Ca*° no contribution from inelastic scattering is expected and the nucleus is so 
heavy that the recoil loss was not large enough to separate the scattering from diffe- 


rent isotopes. 


or 
Or 
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3. Results es ce : tena 


In Tables 2-9 are given the differential cross sections, the Boats angles, the 
rms angular resolution and the mean proton energy, all in the iaboratory system. 
The quoted errors in the cross sections are combined from the following sources: e 

~ (1) Counting statistics. 

(2) Relative errors in the graphical method of determining areas. 

(3) Relative errors in the calibration curve. 

(4) Errors due to the subtraction of inelastically scattered protons. 

(5) Uncertainty of the relative spectrometer field measurement. Relative errors 


due to an incorrect angular dependence of the p—p differential cross sections of Table 1 | 


were not considered. 
In addition there is an uncertainty of the scattering angle of + 0.05° due to varia-_ 
tions in the incoming beam direction and difficulty of positioning the spectrometer at 
the correct angle. 
The absolute errors are caused by the uncertainty of the p—p cross section data 
given in Table 1, and should be less than about 10%. 


Table 2. Lithium. 


Rms angular Differential Mean proton 
resolution cross section energy 
A 6,qy (degrees) | ©),, (mb/sterad) | E,,, (MeV) 


Scattering angle 
01,» (degrees) 


5.34 191.6 +64 
7.84 204.7 3.2 
10.34 146.6 2.7 
12.82 0.52 93.3 2.3 181.9 
15.32 55.4 Ll 
17.80 31.27 0.57 
20.32 17.02 0.37 
25.32 5.866 0.097 
27.82 3.455 0.069 
30.31 2.983 0.041 
32.83 1.623. 0.081 
35.81 1.099 0.029 
7. 0.730 0.038 
40.31 ee 0.482 0.012 182.6 
45.30 0.2069 0.0056 
50.31 0.0728 0.0031 
55.31 0.0303 0.0020 
60.31 0.0108 0.0011 
65.31 0.00587 + 0.00073 


0.00470 0.00041 |} 
0.00519 0.00040 |) 
0,00276 0.00030 |} 
0.00105 + 0,00016 


ane (i 
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“£ * 4, Optical model analysis 
4.1. General 
i A two-body local optical model interaction of the form 
E- V (1) = Vout (1) + U fy (0) +4 Wfachy (0) + (1-2) 9 (09) 
ee 2 
q (EY [othe sew th] 2-0 a) 


Mn dr dr |r 


was used where the form factors are given by 


Pe 
+ - 


a EN 
g (r) =exp [-(" 2) | and R,=r,N1'. 


a 


-N denotes the number of nucleons in the target nucleus. Voou (7) is the Coulomb 
potential from a uniformly charged sphere of radius 1.3 -N7!3- 108 cm, U and W are 
the real and imaginary parts of the central potential and Us and W, the corresponding 
spin-orbit potentials. The imaginary or absorbing part of the central potential has 
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; 2.609 0,083 
28.3 - 204 0.035 
30.27 I~ o.59 1.459 0.044 
32.32 “ 0.902 0.018 
32.81 0.835  0.030_—‘| | 
35.29 | 0.5462 0.0077 || 

37.31 i 0.4506 0.0072 
37.80 0.433 0.016 
39.23) FH . 0.3701 0.0056 | 
40.27 0.349 0,012 
42.78 0.283 0.012 _ 
43.31 0.2879 0.0058 | 
45.28 — 0.2472 0.0036 
46.50 | 0.2056 0.0049 
47.31 0.1885 0.0041 
47.76 0.1797 0.0090 | 
48.52 0.1639 0.0046 
50.28 0.1091 0.0024 
51.50 « 0.0961 0.0032 
52.80 0.0718 0.0038 
55.26 0.0440 0.0026 
55.26 0.0437 0.0017 
60.27 0.01448 0.00064 
65.27 aes 0.00513 0.00027 
70.27 . 0.00285 0.00024 
75.27 0.00183 0.00019 
80.28 : 0.00104+ 0.00011 


0.088 


a form factor that is an arbitrary mixture of Saxon-Woods and Gaussian forms, con- 
trolled by a mixing parameter «. The other symbols have their usual meanings. If 
this potential is inserted in the appropriate Schrédinger wave equation the corre- 
sponding observables can be calculated by the methods described elsewhere [20]. 
Except for a few preliminary runs all the calculations in the present work were done 


on the Oxford Mercury electronic computer. 
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Table 5. Aluminium. 
gre eld ee Se 
1 
Scattering angle Rms angular Differential Mean proton 


6 eessex resolution cross section energy 
tan, (Grerces) A 6,,y (degrees) | 0;,,, (mb/sterad) E,,, (MeV) 


:- 7.81 977.2 +195 
10.31 575.6 9.8 
~~ 49,81 257.3 4.4 
15.30 94.18 1.42 
17.80 30.29 0.39 
19.81 16.15 0.22 
20.27 13.27 0.46 
21.31 11.82 0.23 
21.52 11.61 0.20 
22.76 9.69 0.24 
| 23.30 10.54 0.21 
28.70 10.09 0.11 
10.02 0.18 
25.27 0.55 8.72 0.13 0 
26.54 7.515 0,098 
27.79 5.921 0.089 
30.26 3.944 0,055 
32.77 2.302 0.046 
35.25 1.402 0.051 
36.32 1.083 0.021 
37.77 0.818 0.017 
40.23 0.418 0.013 
49.75 0.2286 0.0066 
- 45.29 0.1310 0.0030 
47.76 0.0881 0.0026 
50.25 0,0626 0.0018 
50.27 0.0726 0.0023 
55.31 0.0232 0.0019 
60.29 0.91 0.00886 0,00054|$ 177.7 
65.27 0.00337 0.00056 
70.27 0.00158+ 0.00022 


4.2. Phenomenological approach 


The general form of the potential (1) contains 15 adjustable parameters whereas 
the most commonly used optical potential with all /;’s identical and « = 1 has only 
six. Although the present data on differential cross sections out to high angles, to- 
gether with existing polarization measurements [6, 9-13] and reaction cross sections 
[8], provide more complete sets of experimental data than have hitherto been 
available they would probably not be sufficient for unique values of all 15 parameters 
to be determined in a free search. Therefore, in the phenomenological analysis de- 
scribed below, essentially seven parameters were varied, viz. U, W, Us, Wtf =T., 
15 =1y, @ = M, =A, =A, = a,. A few trials were made with «+ 1, but otherwise the ab- 
sorbing part was always kept in the Saxon-Woods form. 

In previous optical model analyses of proton scattering at high energies using six 
parameters quantitative fits were obtained to the small-angle data [4, 6, 21-23]. If 
the same parameters are used to calculate the differential cross sections corresponding 
to the present data it is found that in all cases they show marked oscillations, whereas 
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21.56 25.07 0.48 . 
22.82 Oia 20.78 0.37 183.6 
24.06 16.69 0.53 

25.32 11.84 0.19 

27.81 5.92 0.11 

30.31 2.848 0.097 

32.80 1.412 0.034 

35.30 a 0.858 0.018 

37.81 0.588 0.017 

40.33 0.438 0.016 

42.83 0.297 0.013 

45.31 0.1942 0.0074 

45.30 0.198 0.018 

47.82 0.1031 0.0084 

50.31 0.0475 0.0022 

52.81 0.98 0.0220 0.0014 179.5 
55.31 0.01314 0.00089 

57.81 0.0127 0.0012 

60.31 0.0102 + 0.0013 


one of the most noticeable features of the present data is the monotonic fall of the 
cross section with angle, apart from a single oscillation around 20-30°. This discre- 
pancy cannot be ascribed to experimental angular smearing or to contamination 
from inelastic processes as particular care was taken in the present work to achieve 
high angular and energy resolution. f 
The first stage in the optical model analysis of the new data is to try to acpi 
the smooth decrease in the differential cross sections. The data for aluminium were — 
chosen for detailed analysis because it is a fairly light nucleus, so the computing 
time is not too long, and because a marked discrepancy between the observed and 
calculated polarizations was found in previous work on this nucleus [6]. ; 
A series of calculations was made by varying the parameters U, W, a,, a, and dy = yl 
in turn, and in each case the marked oscillations persisted. The next variation studied 
was the effect of different radii for the real and imaginary parts of the potential. 
Keeping the radius of the real part constant at r, = 1.25 f (1 f = 10-18 cm), the radius 
of the imaginary part was progressively increased until at r,= 1.75 f it was found 
that the cross section fell smoothly from 25° to 55°. The corresponding reaction cross 
section, however, was too high, but it was found possible to reduce it to the experi- 
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Table 7. Iron. 
ee ee ee ee eee ee ee 


Scattering angle Rms angular Differential Mean proton 
Dia» (degrees) resolution -eross section energy 
& A Oa» (degrees) | 0,,, (mb/sterad) Ey,» (MeV) 
a ee ee | eae eee 


0.0537 0.0040 
0.0279 0.0019 F 
= 0.01028  0.o0097|f 1/55 


0.00251 0.00043 


10.34 ) 889 +22 
12.79 214.1 3.4 
14.14 66.4 cs 

; 15.32 29.30 0.91 

b 15.93 27.73 0.78 
: yee, 26.59 0.69 

27.33 0.71 

17.82 eee 36.36 0.73 
19.01 42.65 0.81 

: 20.33 44.03 0.42 

: 21.51 37.53 0.75 
22.83 25.09 0.31 
24.04 18.78 0.43 
25.33 12.30 0.15 Hoe 
27.80 4.78 0.12 
30.31 3.113 0.057 
32.81 2.042 0.037 
35.31 1.340 0.026 
37.82 0.887 0.020 
40.31 eee 0.555 0.013 
42.81 at 0.2959 0.0101 
45.31 0.1427 0.0068 
47.81 0.0601 0.0040 

- 50.31 0.0417 0.0030 
52.81 0.0331 0.0024 
55.30 J 0.0207 0.0016 


~ 


mental value by reducing both 7, and 7, in the same ratio without losing the qualita- 
tive agreement with the differential cross section. Subsequent alteration of the absorb- 
ing potential W enabled the oscillation around 23° to be fitted to the experimental 
data without spoiling the agreement with the measured reaction cross section. A 
preliminary account of this study of the effects of different radii has already been 
published [24], together with some of the angular distributions obtained. In all this 
work the radii r, and r, of the spin-orbit potentials were kept equal to the correspond- 
ing radii r, and 7, of the central potentials. 

The primary aim of this work was to fit the differential cross section, but it was 
found that at the same time the fit to the polarization was substantially improved. 
In particular, the minimum around 25° is given by the model, whereas previously 
there was only a slight oscillation in the polarization at this angle. The present data 
for aluminium may therefore be qualitatively fitted by an optical potential of the 
form (1) with the imaginary part of the potential extending some 30% beyond the 


real part. )' . 
A potential of this type was then used to calculate the differential cross sections, 
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52.81 0.0338 0.0021 
55.31 0.0169+ 0.0025 | 


polarizations and reaction cross sections for elastic scattering from the other nuclei 
studied in the present work. In each case overall qualitative agreement was found, 
with the observed oscillation around 25° followed by a steady decrease of the dif- 
ferential cross section well reproduced by the model. 
Further adjustment of the parameters enabled the agreement with the experi. 
mental data to be improved for each individual nucleus, and the best fits obtained 
are shown in Figs. 6-21, and the corresponding parameters in Table 10. As the cal- 
culations are rather lengthy (those for gold, for example, had to take into account 35 
partial waves) it was possible to make only about 8-10 calculations for each nucleus. 
The fits shown could therefore certainly be improved if more detailed parameter 
variations were made. Table 10 also shows the values of A, i.e. the measures of the 


goodness of fit, defined by 
E-T\? . 
0-3 (GE) an 
4 


where 7 is the theoretical value and E +0 the experimental, and the summation 
extends over all the experimental points. The values of A obtained are high, but this 
is most likely more attributable to the incompleteness of the present analysis than. 
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Table 9. Gold. 


Scattering angle Rms angular Differential Mean proton | 
Dray (degrees) resolution cross section energy 
A 6,» (degrees) | ,,, (mb/sterad) | E,,, (MeV) 


ee re | 


9.33 1584 any) 
11.57 0.66 400.0 7.6 181.8 
14.07 566 10 
10.32 557 16 
12.81 531 UT 
15.31 293.4 5.6 
c 17.82 59.0 1.3 
19.31 29.96 0.72 
rie 40.8 1.0 
; 21.59 51.9 1.2 
22.81 oy 45.8 Ll 181.5 
25.31 19.21 0.44 
27.81 6.57 0.22 
29.07 6.55 0.21 
30.34 6.10 0.22 
31.57 a. 27 0.18 
32.83 J 3.84 0.12 
30.32 5.86 0.10 
35.31 1.952 0.041 
37.81 1.056 0.027 
40.31 0.661 0.020 
:. 42.80 1.29 0.412 0.015 182.7 
i 44.05 0.321 0.013 
45.31 0.260 0.012 
47.80 0.1209 0.0079 
50.31 0.0842+ 0.0064 |. 


to any inadequacy of the seven-parameter model. In calculating A no allowance was 
made for the effect of the uncertainty in angle and the angular resolution of the 
detecting apparatus. If this is included A is substantially reduced (cf. Figs 20 and 21 
where the effect of the angular resolution is shown). New computer programmes are 
at present being written, and will be applied to improve the optical model analyses 
made here. 

To illustrate the sensitivity of the distributions to individual variations of the 
optical model parameters, a series of calculations was made with different values of 
U, W, Us, Ws, 71, 72 and a for aluminium. These are shown in Figs. 22-35. It is seen 
that the differential cross sections are mainly affected in the following way: U and Us 
displace the curves vertically without altering their shapes, a determines the average 
slope and the other parameters regulate the diffraction pattern. It is also interesting 
to note the effects of parameter variations on the polarization. The radial para- 
meters /,, 7, and a leave the first polarization maximum essentially unaffected but 
change the slopes of the polarization oscillations in a characteristic way. The potential 
depth parameters on the other hand alter the first polarization maximum but leave 
the oscillations unaffected. The values of the reaction cross sections are also given in 
Figs. 22-35. As might be expected the reaction cross section is very insensitive to the 
parameters which determine the real parts of the potential. 


38:6 559 


10° hae 


Fig. 6. Differential cross sections in CM system for elastic scattering from lithium. All det 
about the theoretical curves are given in Tables 10 and 11. , 


ee EST. hit 
_s==-- Calc. parameters. 


~05 Exp. at 181 MeV [13]. 
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10° 20° 30° 40° 


0° 50° 60° 70° M cm 
Fig. 7. Polarization of protons elastically scattered from lithium. All details about the theoretical 
curves are given in Tables 10 and 11. The energy difference between experimental points and 
theoretical curves has in all polarization diagrams been approximately corrected for by trans- 


forming the angular scale for the experimental points in such a way that ? VB was kept constant 
for each point. 
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| eross sections in CM system for elastic scattering from beryllium. 


Best fit. 
—---—- Calc. parameters. 
} Exp. at 161 MeV [13]. 
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Fig. 10. Differential cross sections in CM system for elastic Scnitecing fronicaereaee 
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Best fit. 
—---— Calc. parameters. . 
$ Exp.at 173 MeV [fi]. 
$ Exp. at 155 Mev [9]. 
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Fig. 11, Polarization of protons elastically scattered from carbon. 


562 


Al 


Best fit. 
-—---Calc. par,volume abs. 
te Calc. par, surface abs. 
$ Exp. at 175 MeV, [12] and unpublished. 
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Fig. 13. Polarization of protons elastically scattered from aluminium. — 
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Fig. 14. Differential cross sections in CM system for elastic scattering from calcium. ; 
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Best fit. 
—-—-- Cale. parameters. 
$ Exp. at 175 MeV [ii] 
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Fig. 15. Polarization of protons elastically scattered from calcium. 
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‘Fig. 16. Differential cross sections in CM system for elastic scattering from iron. 
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Best fit. 
——--= Calc. parameters. 


} Exp. at 179 MeV [13]. 
$ Exp. at 155 MeV [10] . 
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Fig. 17. Polarization of protons elastically scattered from iron. 
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. 18. Differential cross sections in CM system for elastic scattering from indium. 


Best fit. 
-—---— Calc. parameters. 
$ Exp. at 160 Mev [6] 
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Fig. 19. Polarization of protons elastically scattered from indium. 
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Best fit. 
Deena Best fit angular res. folded. in. 
-—---Calc, parameters. 
$ Exp. at 160 MeV [6] 
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Fig. 21. Polarization of protons elastically scattered from gold. 
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Fig. 22. The effect of varying potential parameters on the differential cross sections for alumin 
The parameters of the solid curve is: U = —16 MeV, W = —8 MeV, Ug=2.5 MeV, Wg= —1.0 Me 
7, =7, =1.15 f, rp =r, =1.50 f and a =a, =a, =a = Ay = 0.55 FL i¢ 
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Fig. 23. The effect of varying potential parameters on the polarization for aluminium. eae 
of the solid curve are given in the text to Fig. 22. 
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Fig. 25. See text to Fig. 23. 
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Fig. 26. See text to Fig. 22. 
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Fig. 27. See text to Fig. 23. 
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Fig. 28, See text to Fig. 22. 
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Fig. 29. See text to Fig. 23. 
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30° eine 
Fig. 30. See text to Fig. 22. 
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Fig. 31. See text to Fig. 23. 
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’ Fig. 32. See text to Fig. 22. 
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Fig. 33. See text to Fig. 23. 
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Fig. 34. See text to Fig. 22. 
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Fig. 35. See text to Fig. 23. 
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4.3. Calculation of potential parameters 


_ The six-parameter local optical potential was first expressed in terms of the for- 
ward amplitudes of the nucleon-nucleon scattering matrix and the nuclear density 
by Riesenfeld and Watson [25]. It was subsequently shown by Bethe [26] and by 
McManus and Thaler [27] that the inclusion of the angular dependence of these ampli- 
tudes will imply different radii for the potential terms. Extensive studies along these _ 

_ lines have been made by Cromer [28]. In the present parameter calculation the method 
described in detail by Kerman, McManus and Thaler [7] was followed. 
- The nucleon-nucleon scattering matrix was assumed to have the usual form 


M (q) =A (q) + B(Q) Gyn Gon + € (Q) (G1n + Gon) + E (q) 649 Gag + F (9) Orn Gap, (2) 


CE 


where fig is the momentum transfer and the amplitudes depend on isotopic spin 
_ which will be indicated by indices. For even-even nuclei the B, H and F terms cancel 
exactly when the spin average over the target nucleons is taken; they were always 
_ neglected in the calculations below. The exact number of neutrons and protons in the 
- nuclei was taken into account using the following spin averages for A and C: 


_ N-22 
A (q) =2 (49 +34) + Gay (Ay— Aj), (3) 
aia 22 
Cl) = 2 (Op +80,) + ~~ (CoO) (4 


> 


where N is the number of nucleons and Z the charge of the target nucleus. Instead of 
C(q) it is more convenient to use O(q) given by 
C(q)-k sin 6 = ky O(a), (5) 
where 6 = the scattering angle in the many-body CM system, 


ii =the momentum of the incident proton in the many-body CM system, 
hk, =the momentum of the incident proton in the two-body CM system. 
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: wher pa Wiie the nucleon mass and ar @) the ett form oe aati 


a be ee 
F(q)=Je*"o(r)dr. 


The differential cross sections obtained after solving the ee. equation 1 u 
the potential (6) are to be multiplied by (N’/LN — 1])? before being compared 
experiment. The calculation of the potential parameters is illustrated below for the 
real central term. (The other parameters were obtained in a corresponding manner.) _ 
- There are three parameters to be determined: U, R, and a,. The following equations” 
are valid to a very good approximation: 


[Of (dr= "ZR, (23+ aPal) -U, (10) 
[Uf (r) 2dr =(3Ri+ ga ai) - SUF, (r)dr, (11) 
JU fp, (r)rtdr = (8 Ri +8 Rina} +8 atat)- SUP (r)dr. (12) 


Integration of the real part of equation (7) gives 


Pret’. tiidr = —4a-" (W-1)Re{4(O}, (13) 


P'"(0) __, Re {A" (0)} 
F(0) ~ Re{A(0)} 


PRe{Ve(n}redr=|—3 ]-fRe(reimhar, (14) 


FY (0), F’’ (0) Re{A” (0)} 

Re{V,(r)} rdr=|5--——O) + Re{A" (0) _ 

Base | F(0) *’ FO) Reap} 
5 Re {4% (0) 


is Soe [Re{V,(n}dr. Cal 


Setting (10) = (13), (11) = (14), and (12) = (15) gives three equations from which U, 
R, and a, can be determined if F'(q) and A(q) are known in a series expansion about — 
q = 0 up to terms of the fourth power. From equation (9) F(0) =1, and the deriva- 
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Table 11. Calculated optical potential parameters. 


All lengths are given in fermis and all energies in MeV. 


Al volume} Al surface 


ae absorption jabsorption oo He | i as 
ee lie fs ee | eee ee eee 
PT, 1.609 1.509 1.418 1.241 1.241 1.202 1.163 1.115 1.109 
a, 0.465 0.485 0.509 0.618 0.618 0.688 0.690 0.710 0.746 
U |— 8.45 |—10.39 |—12.88 | —19.44 —19.44 |—21.55 |—24.63 |—29.51 |}—81.08 
Te 1.475 1.390 1.310 1.184 0.994 1.153 1.134 1.103 1.103 
| ay 0.450 0.471 0.496 0.597 0.949 0.665 0.662 0.673 0.704 
W+\—10.44 |—12.52 |—14:82 | —20.91 = 23-060) 922716) 25.33 |— 30.01 | 31.39 
om 1.000 1.000 1.000 1.000 0.198 1.000 1.000 1.000 1.000 
r2 0.946 
b : 1.135 
Eats 1.299 1,245 1.198 1.125 1.125 1.120 1.105 1.091 1.098 
as 0.437 0.454 0.477 0.570 0.570 0.637 0.627 0.624 0.647 
Us; 1.98 2.33 2.50 3.02 3.02 3.05 3.34 3.70 3.73 
Te 1.306 1.252 1.205 1.130 1.130 1.125 1.108 1.092 1.098 
CA 0.466 0.482 0.503 0.591 0.591 0.655 0.647 0.647 0.670 


Wa 0.54 {/— 0.60 |— 0.61 = OORT SOT) eee Os 0) Na OLS Om In a1 OGM ieemal lO) 


_ tives of F'(q) are closely related to the second and the fourth moments of the nuclear 
density. These moments were found from electron scattering experiments [14]. 
Interpolation was used for nuclei for which no electron data are available. The values 

- for A(q) were obtained from the Gammel-Thaler nucleon-nucleon potentials by inter- 

_ polating to 180 MeV between the values tabulated by Kerman et al. [7]. The para- 

* meters obtained in this manner are given in Table 11. They have been used to cal- 
culate the differential cross sections and the polarizations, which are shown in Figs. 6— 
21 together with the best fits from the phenomenological parameter search. In the 
ease of aluminium two kinds of absorbing potentials were used, one with a pure 
Saxon—Woods form factor and the other with such a combination of volume and 
surface absorption that the volume and the second and fourth moments of the poten- 
tial were still given by equations (13)-(15); i.e. judged by the method used for deriving 
potentials there should be no difference between the two types of imaginary poten- 
tials. 


4.4. Discussion 


The overall agreement between experimental results and the best fits of the pheno- 
menological potential suggests that the optical model is applicable to the large-angle 
scattering of high-energy protons by nuclei, but this agreement can only be obtained 
if the radii of tle potential terms are allowed to be different. It is seen from Table 10 
that r,~1.5 r, for the lightest nuclei and r, + 1.2 7, for the heavier ones. This con- 
clusion may be modified when more data are available or when more extensive 
parameter searches are carried out. In particular, when 7; and 7, are allowed to differ 
from 1, and r, respectively, one might find best fits with parameters that are substan- 
tially different from those given in Table 10. Such investigations will, however, have 
to await the completion of faster computer programmes. 

The scattering predicted from the potentials with calculated parameters also shows 
an overall agreement with experimental data. However, the oscillations in the differ- 
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nonlocal and thus it does not necessarily imply that the local approximation is 
invalid. The reaction cross sections given by the potentials of Table 11 are compared 
to experimental data in the preceding paper [8]. From this comparison one may 
conclude that the magnitude of the imaginary potential W should be reduced by 
about 20 % for a reasonable agreement between theoretical and experimental reaction 
cross sections. Such a change would also tend to improve the fits to the differential 
cross sections at small angles. Generally it would be interesting to see if slight devia- 
tions of the potential parameters from the values given in Table 11 are able to produce 
an optical potential that predicts a scattering which more closely agrees with experi- 
mental results. 


ACKNOWLEDGEMENTS 


We want to thank Professor The Svedberg for his interest and support, Mr. P. Isacsson for 
the design of the fast electronic circuits, and Mr. A Svanheden and the cyclotron crew for valuable 
help during runs. The experimental part of this work was supported by the Swedish Atomic 
Research Council. 

We are grateful to the Director of the Computing Laboratory, Oxford, for the use of the com- 
puter, Dr. B. Buck and Mr. R. N. Maddison for their work on the computer programme, and the 
technical staff for their assistance. Thanks are-also due to Dr. K. Parker for performing some 
preliminary calculations on the A.W.R.E. IBM 709 computer at Aldermaston. 


The Gustaf Werner Institute for Nuclear Chemistry, University of Uppsala, Sweden, and the 
Clarendon Laboratory, Oxford, England. 


REFERENCES 


Taytor, A. E., Reports on Progress in Physics, Vol. XX, p. 86, 1957. 

Farssner, H., Ergebnisse der exakten Naturwissenschaften, Vol. 32, p. 180, 1959. 
Squires, E. J., Progress in Nuclear Physics, Vol. 8, p. 47, 1960. 

Bsorxuunp, D. F., Buanprorp, I., and Fernsacn, S., Phys. Rev. 108, 795 (1957). 
Witson, R., Phys. Rev. 114, 260 (1959). 

JoHansson, A., TrBELL, G., ParKER, K., and Hopason, P. E., Nucl. Phys. 21, 383 (1960). 
Kerman, A. K., McManus, H., and THater, R. M., Ann. Phys. 8, 551 (1959). 

Jonansson, A., SVANBERG, U., and SunpBERG, O., Arkiv Fysik 19, 527 (1961). 

ALPHONCE, R., JoHANSSON, A., and TIBELL, G., Nucl. Phys. 4, 672 (1957). 

ALPHONCE, R., JoHansson, A., and Trpety, G., Nucl. Phys. 4, 643 (1957). 


SS eA SN a OSS 


_ 


578 


rT nee 


Tryckt den 26 april 1961 


Uppsala 1961. Almqvist & Wiksells Boktryckeri AB 


579 


Qa" > a the is 
Acid y Sa Toth aan 


4 


ae 
ae 


